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Abstract 



Let be a two-dimensional heat conduction body. We consider the problem of 
determining the heat source F{x,t) = ip{t)f{x,y) with ip be given inexactly and / 
be unknown. The problem is nonlinear and ill-posed. By a specific form of Fourier 
transforms, we shall show that the heat source is determined uniquely by the minimum 
^ ■ boundary condition and the temperature distribution in fi at the initial time t = and 

CN ! at the final time t = 1. Using the methods of Tikhonov's regularization and truncated 

integration, we construct the regularized solutions. 
MSC 2000: 35K05, 42B10, 65M32. 
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1. Introduction and main results 

^- 

j^ ■ Let fi = (0,1) X (0,1) be a heat conduction body and u{x,y,t) be the temperature 

in Q. We consider the problem of determining a pair of functions {u, f) satisfying the system 

ut- Au = (p{t)f{x, y), te (0, 1); (x, y) E n, 

M^(0, y, t) = M^.(l, y, t) = Uy{x, 0, t) = Uy{x, 1, t) = 0, (1) 

uil,y,t) = 0, 

subject to the initial datum and the final datum 

uix,y,0) = go,u{x,y,l) = gi. 

Here, ip G L^(0,T) and go, gi G L'^{^) are given inexactly. 

This is a case of the problem of finding heat source -F(^, t, u) satisfying the heat equation 

Ut — Au = F, 
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for ^ is the spacial variable. The problem has been investigated intensively for the last three 
decades by many authors. Because the problem is severely ill-posed and difficult, many 
preassumptions on the form of the heat source are required. Recently, in [U El El IH E] , the 
authors reduced the heat source F to the function that has separated form 

F(e,t,w) = y.(t)/(e) 

with ip is unknown. Then, the authors in [21 [7j studied the problem in the case that the 
heat source F{^,t,u) = (f{t)f{^) with / is unknown. From some assumptions of <p, the 
author used the Fourier transform and truncated integration to regularize the problem with 
nonsmooth data. However, in [7], the Cauchy datum u{x,y,t) is given on two parts of the 
boundary, i.e the boundary x = 1 and the boundary y = 1, say. In the present paper, the 
Cauchy datum is given only on the boundary x = 1. The requirement of the Cauchy datum 
is minimum because if the condition u{l, y, t) is omitted then the uniqueness of the solution 
of the problem cannot hold. For example, we consider the system 

ut- Au = ip{t)f{x,y), 

M^(0, y, t) = Ucc{l, y, t) = Uy{x, 0, t) = Uy{x, 1, t) = 0, 

subject to 

u{x,y,0) = go,u{x,y,l) = gi. 

This system has not the uniqueness property. Indeed, if 

do = 9i = 0, ip(t) = TTCos{nt) + 27r^ sin(7rt) 

then this system has, beside the trivial solution [u, /) = (0, 0), a nontrivial solution 

u{x, y, t) = sin(7rt) cos(7rx) cos(7r|/), /(x, y) = cos(7rx) cos(7r|/). 

In fact, by a specific form of Fourier transforms, we shall get 

Theorem 1. (Umqueness) Let ui,U2 G C\[0,l];L\n)) f] L'^{0,1] H^{^)), fiji e L'^iP)- 
If {ui, fi) (i = l,2j satisfy the system (1) with go,gi G L'^{Q) and ip G L^{0, 1)\{0} then 

(uiji) = (m2,/2)- 

We also have a regularization result. Using the Tikhonov regularization and truncated 
integration, we can construct a regularized solution for all (p ^ 0. 

Theorem 2. (Regularization 1) Let {uexjex) G {C^{[Q,l]]L'^{yt)) f\ L'^{{),l]H'^{VL)),L'^{yL)) 
he the exact solution of the system (1) corresponding the exact data goex,giex G L'^{Q) and 
ipex G L^{0, 1)\{0}. Let go£,gie G L'^{fi) and (p^ G 1/^(0, 1) be measured data satisfying 

llfi'Oe ~ fi'0ex||/,2(Q) < S, \\gie — giex\\L'^(^Q) ^ ^, W'-Pe ~ V^exH ^1(0,1) ^ ^• 

From {goe,gie,^e}, we can construct a regularized solution fis G L'^{R^) such that 



Moreover, if f^x G H^{Vt) then for each j3 G (0, 1/2), there exists Eq > (depended on 
Uex, ^ex dnd P) such that 

Wfie - fexWl^n) < 4vr-' (l + 9 WfexW^.^nyj M^^'T^ 

for < e < Eq. 

A smoother regularized solution will be given if ip satisfies the following condition 

(H) There exist A G (0, 1) and Cq > such that either ip{t) > Co for a.e t G (A, 1) or 
(f{t) < -Co for a.e t G (A, 1). 

We note that (H) will be satisfied if ip is continuous at t = 1 and v'(l) 7^ 0. Under (H), 
one has 

Theorem 3. (Regularization 2) Let (uexjex) e {C^{[0,1]; L'^i^)) ^ L^{0,1; H^{n)), L'^{n)) 
be the exact solution of the system (1) corresponding the exact data goex,giex G L^(f2) andipex 
satisfying (H). Let ge = {gos^Qie) ^ (-^^(^))^ '^'^'^ Ve ^ -^"^(0, 1) he measured data satisfying 

llfi'Oe ^ 9Qex\\L'^{yi) ^ ^j Hfi'le — 9lex\\L'^{yi) ^ ^j We — V'ex |Ili(o,1) ^ ^• 

From {gQe^gieWe}, we can construct a regularized solution /2e G L'^{R^) such that 

\mi\\f2e- fexWLVn) = 0- 



e^O 



Moreover, if f^x G H^ {Vt) then there exist 7 > (depended on ipex) and £0 > (depended 
on Uex, Vex) such that 

11/2. - /e.||i.(f,) < 47r-l (||/e.||ix(^) e^ + 9 ||/e.||HX(^) S^") 

for < £ < eo • 

The remainder of the paper is divided into three sections. In Section 2, we shall give 
some notations and preparation results. The main results will be proven in Section 3. In 
Section 4, a numerical experiment will be given to illustrate our approximation. 

2. Notations and preparation results 

First, we have 

Lemma 1. If u e C^{[Q,1]]L'^{Vl)) n L'^{Q,1]H'^{Vl)) and f G ^^(fi) satisfy the system (1) 
then for all {a,n) E R x Z we have 

g{x,y) — e^^°' "^" '^ 'g{x,y)j cos{ax) cos{mTy)dxdy 

1 

= f e("'+"'"')(*-^VW^^- / fix, y) cos(aa;) cos{mxy)dxdy. 
n 
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Proof. Getting the inner product (in L^{Q)) the first equation of the system (1) with W = 
cos{ax)cos{nny) , we have 

— / uWdxdy + (a^ + n^vr^) / uWdxdy = (pit) / fWdxdy. 
a. n n 

We multiply the latter equation with e*^" "''"' ^ )(*~i) to get 

g(a2+„2^2)(,_i) f^^^^^y] ^ g(a2+n2.2)(i-l)^^^^ f fWdxdy. (2) 



d 

di 



n / n 



Integrating (from to 1) the latter equality with respect to t, we shall get the desired 
result. D 

Definition 1. For all ip E L^{0, 1), we set D{ip) : R x Z ^ R 

1 

D{ip){a,n) = /"e("'+"'"')(*-iVWt^^, 



and 

B{ip,r, 0") = {« G {—r,r) \3n E Z H {—r,r), \D{ip){a,n)\ < a} . 

Let A be a subset of R. From now on we denote by m{A) the Lebesgue measure of A. 
Using the idea in [S] (see Theorem 4), we have the following result. 

Lemma 2. Let ip e ^^(0, 1), p G (0, 1/2) and q> qi>0. 

[i) Assume that ip ^ 0. Then D{ip){a,n) ^ for all n E Z and for a.e a & R. Moreover, 

there exists ^o > (depended on ip, q and [3) such that 

miB{ip,Ru,e''))<R^} 

for < e < Eq and Ri^ := (ln(5"^))^. 

[ii) Assume that ip satisfies (H). Then there exist 7 > (depended on ip) and Eq > Q 

(depended on ip, q and qi) such that 

m{B{ip,R2e,e''))<e^ 
for < e < Eo and i?2e := e^^^^^. 
Proof Because (p G L^(0, 1)\{0}, the map (p : C ^ C 

1 1 

-2 / „2t. 



f e''ip{t)dt = f e-''ip{l-t)dt 



(p[z) = e 



is a nontrivial entire function. Hence, for each n E Z, the function 

1 

0„(z) = ^{z^ + n\^) = /"e(^'+'^'^')(*-iV(t)^i 





is also a nontrivial entire function. 

(i) For each n E Z, since the zeros set of 0„ is either finite or countable, D{ip){a,n) = 

0n(«) 7^ for a.e a E R. Hence D{(f){a, n) ^ for all n E Z and for a.e a E R. 

To estimate the measure of the set B{ip,r,a), we shall use the following result in [9] 
(Theorem 4 of Section 11.3). 

Lemma 3. Let f{z) be a function analytic in the disk {z : \z\ < 2eR}, |/(0)| = 1, and let t] 
he an arbitrary small positive number. Then the estimate 

ln|/(z)| > -ln(^^).ln(M/(2ei?)) 

is valid everywhere in the disk {z : \z\ < R} except a set of disks {Cj) with sum of radii 
'Y^Tj < TjR. Here Mf{r) = max |/(2;)|. 

\z\=r 

Since 0^0, there exists a^ G (—1,0) such that |0(ao)| = Ci > 0. For each n E Z, 
I'^l < Riey we put Zn = iy/u'^n'^ — Oq. Then \zn\ < irRie + 1 and |0n(-2n)| = |0('^o)| = C*i- We 
have that 

n[Z -\- Zn I 



^Jz) :-- 



Ci 



is an entire function and |\1'„(0)| = 1, moreover for all z E C, \z\ < 2eR, n E Z, \n\ < Ri^, 



\^n(z)\ = —. 



^((.+.„)2W^^)(i-l)^(^)^^ 



i 

< ^. /"e(l"l+l^"l)'|(^(t)|rft<e 
<^i J 






For e > small enough (depended on ip, q and (3) and for each n E Z, \n\ < R^, applying 
Lemma 3 to i? = (1 + vr)_Rie + 1 and rj = ^^^^ , we obtain that 



ln|^„(z)| > -ln(l^).ln(M^„(2ei?)) 
V 



> - [\n{R) + 2 ln(i?i,) + ln(75e= 



> -3 ln(i?i,)(2e(l + tt) + tt + lyRi + In f — 

1 



(2ei? + nRis + if + In 



1 



Iv^IIli(o.i ) 



> -gln(£ ^) +ln( — 



'■"c; 



for all \z\ < R except a set of disks {B{znj,rnj)}j^j„ with sum of radii 

J2 ^"i <vR= v^- 



i&Jn 



5R 



Consequently, for e > small enough and for each n ^ Z,\n\ < Ri^ we get 
|D(v9)(«,n)| = |0„(a)| = C,. |vl/„(a - ^„)| > C^.^e'^ = e'^ 
for all \a\ G [—Rie^Ru] except the set U B{znj + Zn,rnj)- So 

B{ip, Ris, e^) C U U {^nj - rnj, inj + Tnj) 

neZ,\n\<Ri^ j£j„ 

with C,nj = Re{znj)- Thus 

for e > small enough (depended on (/?, g and /9). 
{ii) Note that 



L'(v9)(a,n) 



1 







A 



A 



(a^+n^.^)(X-l) \^^t^\dt+ I e^°'"+^"^"^^'-'-^Codt 







y _ (a2+n2^2)(A-l) II II , ^ 



1 _ g(a2+„V2)(A--l) 



Therefore, there exists a constant -Ri > (depended on v^. A) satisfying for either \a\ > Ri 
or \n\ > Ri that 

2(q;^ + n^TT"^) 
Consequently, for e > small enough (depended on ip, q and gi) and for all {a,n) G 

{-R2e,R2e?\i-Rl,Rl?,^eget 



D{^){a,n)> 



Cn 



> 



Cn 



Cn 



-.e"^ > e^. 



2(a2 + n%2) - 2(1 + Tr^)Rl 2(1 + 7r2) 

Now we consider only (a, a) G (— i?i,-Ri)2. Let ao, -Zn, ^n as in (i) and put R2 
[1 + 7r)i?i + 1. Then for all z eC, \z\ < 2ei?2, we get 






^({z+z„f+n^n'^){t-l) 



^{t)dt 



Ci 



where C2 > 1 be a constant independent of n. 

For £ > small enough, applying Lemma 3 to R = R 



'2, ^ 



Qi 



(4i?i + 2)i?2 



and 7 



> 0, we get 



2/n(C2; 

ln|^„(^)| > [7. ln(£) - ln(75i?2e')] . ln(C2) > 27^(^2) ln(£) - In(Ci) = ln(Ci- 
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l^qV 



for all \z\ < R2 except a set of disks {B{znk,i^nk)}keKn with sum of radii 



"^ r-nk < 11R2 



keKn 



4i?i + 2 



Consequently, for e > small enough and for each n & Z,\n\ < Ri, we have 

\D{ip){a,n)\ = |0„(a)| = C,. |^„(a - z„)| > 5^^ > e^ 
for all a G [— -Ri,i?i] except the set U B{znk + Zn^Tnk)- Thus 

B{ip,R2e,e'^) C. U U (C„fc-r„fe,^„fc + r„fc) 

neZ,\n\<Ri keKn 

with ^nk = Re{znk)- Therefore, 



m 



iBiip,R2e, €''))< Yl I] 2r„, < (2i?i + 1).2. 

-Ri<n<Ri keKn 



4i?i + 2 



for £ > small enough (depended on ^9, q and jS). 
The proof of Lemma 2 is completed. 

Definition 2. For eac/i w G -^^(fi), iL'e se^ G{w) defined on R x Z by 

G{w){a,n) = / w{x,y) cos{ax) cos{mry)dxdy. 



n 



Lemma 4. {i) For each v G -^^(0, 1) we have 



v{x) cos{ax)dx 



da = IT \\v\ 



L2(0,l) ■ 



(m) For each w G L'^{Q) we have 



00 
00 „ 

J: J\GH(a.n)fda 



-K \\W\ 



L^{Q) ■ 



Proof, (i) Putting 

(vix), a;G(0,l), 
v{x) = < v{—x), X G (—1,0), 
[0, x^(-l,l), 

then the Fourier transform of v(x) is 



F(f)(Q;) := / v{x)e ^°'^dx = 2 / f (x) cos(aa;)(ia;. 





Using Parseval equality, we get 



v{x) cos{ax)dx 



[a) For each n E Z, we put 



da = -\\F{v)\\l2^j^) = -|klli2(K) =T^\Hl2^o,i)- 



hn{w){x) = I w{x,y)cos{n-Ky)dy. 





Applying (i) to v = hn{w), we get 



\G{w){a,n)\ da 



hn{w){x) cos{ax)dx 



da = n \\h„(w 



n\'^)\\L'2{0,l) ■ 



On the other hand, since w G L'^{Q), for a.e x G (0, 1) we have w{x, .) G L'^{0, 1) and its 
Fourier-cosin series corresponding to variable y is hn{w){x). Using Parseval equality, we get 



[x)\ = WW 



{^^■)\\l^o,i)= \'^i^^y)\ dy, a.e x G (0, 1). 





n=—OD 

Therefore, 

^ / \G{w){a,n)fda = 7r J^ ||^n(w)||5^2(o,i) = tt ^ / |/i„(w)(x)|^(ix 



n \n=-oo / n \ n 



\w{x,y)fdy dx = 7i\\w\\l2(^^y 

V"— ^ / \o / 

The proof is completed. 

To prove the regularization results, we need one more preparation. 
Definition 3. For each w G L'^{Q) and r > 0, we set 

CO 
p OO « 

,(w.r)=j: J \G(^)(a,n)fda+Y. J \GM(a,nf da. 



n 



\n\>r_^ 



\a\>r 



Lemma 5. For each w G L'^{Q), we have lim fi{w,r) = 0. Moreover, if w E H^{Q) then 



n{w,r) < - + 



8 27r 



8 



\w\ 



i/i(f7) • 



Proof. For each w G L'^{Q), applying Lemma 4 we obtain 



E 



G{w){a,n] 



da = 11 \\w\\j^2iQ\ < +c>o. 



It implies that lim fi{w,r) = 0. 

r— »+oo 



Now, we consider w G H^{Q). Since 



w{x, y) cos{nny)dy 



w{x,y) 



sm{miy) 



nj/=i 



nn 



y=0 



dw , , sini nny) 1 f dw 

-^{x,y) dy- 

oy rni 



nil J dy 





{x,y)iAn{miy)dy, 



we get 



G{w){a,n) = / ---{x,y)cos{ax)sm{m[y)dxdy = G{—-){a,n). 

nn J oy rni ay 

n 

Consequently, 

oo oo 

— oo 
_. oo 

2^2 Z^ 



|n|>r_ 



\n\>r 



Gi^-)icy,n) 

oy 



da 



r^TT^ 



G(^)(a,n) 



r^TT 



dw 



dy 



L'HQ) 



Similarly, we have 



Therefore, 



w{x,y) cos{ax)dx 



w{x,y) 



sin[ax) 
a 



x=l 



x=0 



dw , , sinfaa:) , 

^-(a:,!/) dx 

ox a 







1 



smfa) 1 f dw . . . , . , 

w[l,y) / -——[x,y)sin[ax)dx. 

a a J dx 





1 

G(w)(a, n) = / wil, y) cos(n-Ky)dy G I — — ) (a, n). 

a J a \dx J 
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Hence 



oo „ 

J2 / \G{w)fda 



\a\>r 



< 



/2 -^ / pi \ 2 r 

—da. J2 (/ H^,y)cosimry)\dyj +— ^ / 



OX J 



da 



|a|>r 

4 2 2 



(9w 



dx 



L2{n) 



Noting that 



d 

^(l'^)= / ■^{xw{x,y))dx 



w{x,y) + x—{x,y) ] dx, 



we get 



\wil,y)f< 



w{x,y) + x—{x,y) ] dx 



<2j \ \w{x,y)\' + 





l^(^'^) 






and 



Thus 



|u'(l,.)ll^^2.n.^ <2\\wfr.rr..+2 



lL2(o,l) 



lL2(f7) 



(9w 



9a; 



L2(n) 



^ 2 ll'wll^i(f^) 



^ /* ^ 9'7r 

^ y \Giw)f da<- |h||5,i(^) + — 



dw 



dx 



^1 I \ 



L2(f7) 



In summary, we get 



I ^ 27r \ 2 
/i(w,r) < - + ^ \\w\\hi 






(n)- 



The proof is completed. 



D 



3. Proofs of main theorems 



Proof of Theorem 1 



Proof. Put u = ui — U2, f = fi — /2- Then (m, /) satisfies system (1) corresponding to 
9o = 91 = ^- Let D and G be as in Definition 1 and Definition 2. Applying Lemma 1, we 
obtain 

D{ip){a, n).G{f){a, n) = 0, V(a, n) eRx Z. 



10 



For each n E Z, according to Lemma 2, D{ip){a,n) ^ for a.e a E R. It implies 
G{f){a,n) = for a.e a G -R. So / = because of Lemma 4. Hence, equation ([2|) in the 
proof of Lemma 1 becomes 



/ u{x,y,t) cos{ax)cos{nny)dxdy I 



d j ^(„2+„2,2)(,_i) , u{x,y,t)cos{ax)cos{nny)dxdy I = 0. 



Since u{x, y, 0) = 0, the latter equation implies that for all t G (0, 1) 

u{x, y, t) cos(ax) cos{mTy)dxdy = 



Using Lemma 4 again, we obtain m = as desired. D 

Proof of Theorem 2 and Theorem 3 



Proof. The proof is divided into four steps. 

Step 1. Let D and G be as in Definition 1 and Definition 2. For each g = {qq, qi) G {L'^{Q)) 

we put 

From Lemma 1 and Lemma 2, it follows that 

G{fe^){a,n) = — 

for a.e a E R and for all n E Z. From Lemma 4, we have 

fex{x,y) = -. Y] / ''"' ' cos{ax)da cos{mTy). 

We shall construct the regularized solution /^ by the following formula 

fe{x,y) = -. J^ / -jy^j^——^cos{ax)da\cos{miy), 

where R^ and S^ are regularized parameters chosen later. 
It is obvious that /^ G G{R^) and 

G(/.)(a,n) = X {{{-Rs,Rs?) -^^^ 
where x(^) is the characteristic function of the set A, i.e, 

^^(^>(^)-{;; 111 

11 



To prove that fs approximates fex in -^^(^), "we only have to verify that G{fs) approxi- 
mates G(/e^) inL2(i?2). 

oo oo 

step 2. Estimate E I \Gifs) - GiU)\^ da. 



n=— OO— OO 



We first estimate \G{fe) - GiU)\. If (a,n) ^ {-Re^Ref then \G{f,)-G{f,, 
\G{fex)\- Now, we consider (a,n) G {—RejReY- Putting 

^ _ Jii II I \ 

O3 — maX|||(/9ex|lii(o,l) ) \\90ex\\L2(n) ' llfi'lea;|lL2(f^)J , 

by a direct calculation we obtain 

\H{gex)\ < llfl'lex|lLi(f7) + llfi'Oexllii(n) ^ Sos, 

l-f^(fi'e) ~ -f^(5'ex)| < llfi'le " 5'lex |lLi(f7) + Hfi'Oe — QOexWi^^n) < 2e, 

\D{(Pex)\ < |bex|lLl(o,l) < C'S, \D{(Pe) " ^(v^ex)! < ||^e " Vexh^O,!) ^ ^• 



We have 



\Gif,)-GiU\ 



H{g,)D{^,) H{g, 



D\^,)+6, D{^, 



< 



H{ge)D{Ve) H{gex)D{Lpex) 



H{gex)D{(pex) H{ge 



We shall estimate each term of the right-hand side. We have 
Hige)D{ipe) H{gex)D{(^ex 



D'^{(^ex) +5ex D{(pe 



ex I I '-'ex 



6e 



^ \D{^e)D{Vex)\ ■ \D{Vex)H{ge) - D{^e)H{gex)\ ^ 4 \D{^e)H{ge) - D{^,x)H{ge 



< 



\D{^,x)H{ge)-D{^,)H{g,x)\ , \D{^e)H{ge) - D{^,^)H{g,^)\ 



We get 



and similarly, 



Thus, 



\D{^ex)H{ge) - D{^,)H{g,x)\ 

< \D{^ex)\ ■ \H{ge) - H{g,^)\ + \D{if,) - D{if,.,)\ . \H{g, 

< Cs.\H{g, - gex)\ + e. \H{gex)\ , 

\D{^,)H{g,)-D{^,,)H{g,,)\ 

< |Z}(^ex)| . \Hige) - H{g,^)\ + \D{ip,) - D{ip,.,)\ . \H{g, 
+ \{D{^s) - D{^,x)) {H{ge) - H{g,x))\ 

<2C^.\H{g,-g,,)\+e.\H{g,,)\. 



H{g,)D{^,) H{g,^)D{^, 



D^{V,)+5e D^{^,x)+5e 



< 3C3. \H{g, - fc)| (5-1 + 25. \H{g,x)\6;\ 
12 



We next consider 



D\^,,) + 8,- 



If \D{ipex)\ >€" then 

H{g^^)D{ipe^) H{ge^ 



H{ge 



D^ifex) + 4 D{ip, 



D{iPe 



s, 



D^{iPe.)+S, 



■ \Gife. 



< ^r \GUe- 



T2q 



If \D{^^^)\ <e^ then 
H{gex)D{ifex) H{ge 



< x{B{^,^,R„e'^)) .\G{fex)\ < x{B{^,,,R„e'>)) .Wfexhim , 



where B[(pex,Re7^'^) is as in Definition 1. 

In summary, for all (a, n) G i? x Z , we have 

\Gife) - GiU\ < 3C,5;' \H{g, - g,,)\ + 2e5-^ \H{g,,)\ + x{B{^ex.Re.e'^)) WfexW^n) 

+ [x((«,n) i {-Re,Ref) + x{{oc,n) G (-i?„ i?,)^)^^-'^] . \GUex^ 

oo oo 

To estimate ^ / \G{fe) — G{fex)\ c?q;, we square the latter inequality to get 



n=— oo— oo 



\GUe) - G{fex)V < 5QCi6;' \H{g, - g,,)\' + 165^5;^ \H{g,,)\' 

+Ax{B{vex,Re,e'^))\\fex\\\n) 

+ [4x((a,n) i {-R,,R,f)+A5le-"^] . |G(/e.: 

We will consider each term of the right-hand side. Since 

mg^^^t < {\G{gie^x)\ + \G{g,ex)\f < 2 \G{g,,,)\^ + 2 |G(^oe.)|' : 

we get 



oo „ 



Similarly, 



oo °? 



Furthermore, 



oo 
oo „ 

a=— oo 

— oo 

oo 
oo „ 

Y. I X{{a,n) ^ {-R,,R,)^)\G{U)f da = i^{U,R,), 

i=—oo 

— oo 

oo 
oo „ 

J2 J \G{Ufda = 7r\\f, 



lex\\L2(^Q) ■ 
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Thus 



oo „ 

J2 J \G{fs) - G{U)f da 



< lUixCle^f + Q^T^Cle^f + 4m(i?(<^,,, i?„ e^)) \\U£^, 



m 



|2 



Choosing q = 2/5, (5^ = e^^^^, we obtain that 

oo 
oo „ 

||/e-/e.lli.(o)=Vr-^ E / \G{fe) - G{U)f da 

n=—oo ^^ 

— oo 

whereC4 = 208C| + 4||/e.||i2(^). 

Step 3. Proof of Theorem 2. 

We construct /i^ as fe corresponding 

q = 2/5, 4 = 5^/1°, i?, = i?i, = (/n(£-i))^. 
Applying (i) of Lemma 2, for e > small enough (depended on ip^x and P), we have 

Hence, from ([3]), it follows that 

Since limi?i£ = +oo. Lemma 5 implies that lim/i(/e^, i?ie) = 0. Therefore, 

£—+0 e-^0 

lim||/l£-/ex|L2(m =0. 

Moreover, if fex G H^{VL) then Lemma 5 implies 

So we obtain 

Wfie - fexWl.^n) < G.e'l" + 4vr-i ||/e.||L(n) i?r.' + 47r"i (Si^r/ + 27ri?r.') ||/e.||ii(o) 
<47r-i(l + 9||/,,.||i.(^))i?i-' 

for e > small enough (depended on Uex-, ^ex and [3). 

Step 4. Proof of Theorem 3. 

We construct f2e as /^ corresponding 

g = 2/5, 4 = ^'/'°,i?. = i?2e =£'/'. 
14 



(3) 



Applying (ii) of Lemma 2 to gi = 1/3, there exists a constant 7 > such that for £ > 
small enough (depended on ip^x) we get 

Hence, from ([3]), it follows that 

11/. - UA^n) < C^e'/' + Att-'s^ ||/e,.||ii(^) +47r-V(/e.,i?2.). 
Since limi?2e = +00, we use Lemma 5 to get \im fi{f ex, R2e) = 0. It implies 

lim\\f2e- fex\\L2(n) = 0- 
Moreover, if / G H^{Q) then Lemma 5 implies that 

Kfe., R2e) < {SR^I + ^TtR^s) Wfe.Wl^n) = (S^'^' + '^^^'^') \\fe4mm ' 

Thus 

11/. - fe.fmn) < Cas'" + 47r-i WfeAl.^n) ^' + 4vr-^ {%e'l^ + ^i^e'l'') Wf.xfuH^^ 
<47r-i(||/,.||i.(^)e^ + 9||/e..||?,.(^)e^/«). 

for e > small enough (depended on u^x and ^ex)- D 

3. A numerical experiment 

We consider the exact data 

<^ex = TT e , 

g^ex = (cos(7ra;) + 1) cos^ny), 

2 

giex = e~"^ (cos(7rx) + 1) cos{'Ky). 
Then the corresponding exact solution of the system (1) is 

Uex = e"'" *(cos(7rx) + 1) cos(7r|/), 
fex = cos(7rx)cos(7ry). 

For all m = 2,4, 6, 8, ..., we consider the disturbed data 

gom = goex + m~^{cos{m7Tx) - 1) cos(7ry), 
giex = giex + m~^{cos{mnx) - 1) cos(7ry). 

Then the corresponding disturbed solution of the system (1) is 

Um = Uex + m^^e^'' \cos{mnx) — 1) cos(7ry), 

fm = fex + m cos(m7rx) co?,{-ny). 
15 



We get 



I _ II -^ 

\gOm 90ex\\L2(n) ~ 2^' 

I _ II ^""'^ 

\90m 90ex\\L2(^Q-^ - 



Jm Je. 



2m 
m 

L^in) 2 ' 

It means that, when m is large, a small error of data causes a large error of solution. Hence, 
the problem is ill-posed and a regularization is necessary. 

We construct fim and /2m as /i^ and /2e in Theorem 2 corresponding 



e = m 



5. =£9/1° 



m 



-9/10 



V^e — Vm, dOe — dOm, die — Qlmi 



2m — fn 



1/6 



We have 



D{ipm){a,n) 



7rV(l-e-"' 



a" 



if n^ ±1, 

(m — l)7r^e"'^ (1 — e^" ) sin(Q;)(Q;^ + rmr'^) 



2a{a'^ — 7r^)(Q;^ — rn?Ti'^) 
Hence, for each j G {1,2}, we get 



if n = ±l. 



fjm{x,y) = -. Y] 
( «. 



t^j ??i *C ^ \ ty-j m \ __ D . 




H{gm)-D{.Vm) 



\ 



cos(Q;x)(ia 



cos(n7ry) 



2 

TT 

2 

TT 






ij'(fif^)(a, l).£'(v?^)(a, 1^ 
Z}2(^„^)(a,l)+m-9Ao 



\ 



/ 



cos{ax)da 



cos^ny) 



■jm 



J 



[m — l)7r''a;sin(Q;)(Q;^ + rriTr^ 



\ 



v2N-2 



^_ / 2(a2 _ 7r2) (a2 _ ^2^2) ('^4 + ^-9/i0c,4e2u2 (1 _ ^ 
The error estimates between fj^ and fex is given by the following table. 



cos{ax)da 



cos(7ri/). 



/ 
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e = m ^ 


Jim Jex\ L^(n) 


\j2m Jex\ L^(Q,) 


10-2 


0.2499999413 


0.2499999226 


10-6 


0.2495687285 


0.2494516526 


10-12 


0.1181388651 


1.240475046 x 10-2 


10-1^ 


9.038948894 x 10-2 


5.031837329 x 10-^ 


10-20 


5.881288742 x 10-2 


7.385640149 x 10-^ 


10-30 


3.361360762 x 10-2 


1.591595841 x 10-^ 



Table 1. Error estimates between /,>„ and /g. 



Here are some figures of tlie exact solution and the regularized solution with m = 10 
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Figure 1. The exact solution 








Figure 2. The regularized solution /i^ 
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Figure 3. The regularized solution /; 



2m 



Remark 1. fj^ approximates fex well on Q except a neighborhood of the line y = 1. The 
fact is understandable because the value ofu{x, l,t) is not given in the system (1). 
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